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Abstract. Minkowski valuations provide a systematic framework for quantifying different aspects of mor- 
phology. In this paper we apply vector- and tensor- valued Minkowski valuations to neuronal cells from the 
cat's retina in order to describe their morphological structure in a comprehensive way. We introduce the 
framework of Minkowski valuations, discuss their implementation for neuronal cells and show how they 
can discriminate between cells of different types, codes 
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1 Introduction 

Natural phenomena can be understood as causes and con- 
sequences of a continuing interplay between geometry and 
dynamics, or form and function Spatial adjacencies, 

. the specific geometrical features of natural entities, as well 
as the dimensionality of the space where they are embed- 
ded, largely constrain their dynamics and function. For in- 
stance, the proper operation of a mammal's heart depends 
on a suitable diffusion of potentials and waves across the 
heart surface. It is at the central nervous system, however, 
that the interplay between form and function reaches its 
greatest complexity. To begin with, the velocity of signal 
transmission in neuronal fibers (i.e. dendrites and axons) 
depends on the width and length of the fibers. The impor- 
tance of geometry for proper neuronal function is further 
underlined by the fact that neurons are cells specialized 
to establish selective connections. Given the spatial con- 
straints imposed by three-dimensional space, these cells 
have to resort to the most diverse geometries in order to 
implement the required interconnections - as they do in a 
dynamical way during the whole lifetime of an individual. 
As a consequence it is to be expected that morphological 
analyses of neuronal cells provide clues for understanding 
neural dynamics and function. Although a large number of 
investigations have been directed at the neural anatomy 
and geometry (e.g. (HIIHISlEi ) i only the taxonomic orga- 



nization of neuronal cells or the consideration of shape 
abnormalities as subsidies for diagnosis have been con- 
centrated on so far. The study of the shape of neuronal 
cells with objective and mathematically well characterized 
morphometric descriptors is just at the beginning (e.g. [H] 

EKE]). 

In order to be useful tools, such morphological descriptors 
should fulfill the following criteria: First, the extracted 
quantitative features should obey simple transformation 
rules, when the neuronal cell under investigation is sub- 
ject to elementary geometrical transformations such as 
affine or conformal transformations (in- and covariance). 
Second, the obtained measurements should discriminate 
between different classes of neuronal cells. Finally, it is 
important that the estimated features allow for intuitive 
interpretations from the neuroscience point of view. 
Because of their long tradition in modeling and analy- 
sis, mathematics, physics and engineering provide a large 
number of concepts and measures that may be consid- 
ered for studies in neuroscience and neuromorphometry. 
A good example is entropy, which has been used in neuro- 
science because of its close association with the concept of 
information JU|. Other such measures include the fractal 
dimension |TTil2ll3| . lacunarity [TUTTo] . percolation crit- 
ical density |2] and curvature ^5]. Recently, concepts from 
Integral Geometry and in particular the (scalar) Minkowski 
shape functionals were applied in order to characterize the 
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geometry of ganglion cells from the cat's retina |17U18| . 
Minkowski shape hmctionals are particularly interesting 
because they meet the criteria mentioned above: They 
are invariant under rigid body transformations, seem to 
have good discriminative power [TJ\, and can be squared 
with basic concepts from neuroscience. Moreover, they can 
easily be implemented: Usually, the original data are fil- 
tered with methods known from MIA (Morphological Im- 
age Analysis). This preprocessing introduces a free pa- 
rameter, which can later be varied in order to probe the 
morphology at different scales. In previous works, the sin- 
gular points (branching and terminating points) |18| or the 
whole cell outline were dilated, where the dilation ra- 
dius entered as a parameter. For each dilation radius, the 
preprocessed neuron image was decomposed into compo- 
nents (or basic building blocks). The Minkowski function- 
als can then be calculated by counting certain multiplici- 
ties of the basic building blocks. This approach makes use 
of mathematical results from Integral Geometry . 
In this paper, we use extensions of the Minkowski shape 
functionals, viz. the Minkowski valuations, in order to fur- 
ther improve neuromorphometric characterization and ana- 
lysis. These extensions were only very recently investi- 
gated by mathematicians |2()U21U22j and include vector- 
and tensor- valued measures. They are therefore sensitive 
to directional information and also allow for valuable graph- 
ical visualizations. Minkowski valuations have been suc- 
cessfully applied to describe the morphology of galaxies [231 
and galaxy clusters [2"^| . 

In the following, we will illustrate the potential of the 
Minkowski valuations for neuromorphometry by analyz- 
ing a set of ten ganglion cells from the cat's retina. We 
consider two-dimensional projections of the cells. The set 
used .25 includes cells with diverse shapes, correspond- 
ing to a recently revised classification of those types of 
cells |21]. In addition, ganglion cells from the retina exhibit 
branching patterns which are predominantly planar, and 
therefore compatible with the two-dimensional Minkowski 
valuations considered in the present work. 
The article starts by presenting the higher-order Minkows- 
ki functionals and proceeds by illustrating their applica- 
tion to the characterization of neuronal cells. 



2 Minkowski valuations 

Morphometry deals with measures for the content, shape 
and connectivity of spatial patterns ( "bodies" ) . Consider 
a body P in 2-dimensional space such as constituted by 
the pixels of a neuron image (see Figuresn an d|21fbr exam- 
ples). A straightforward way to measure its "content" is 
to calculate its area Vb (P) or - equivalently - to count its 
pixels. The area clearly meets the requirement of motion 
invariance. Furthermore, it is additive; that is, the area of 
the set union P of two bodies P\ , Pi can be decomposed as 
V Q (P) = Vo(Pi) + Vq(P2)-Vo(PiUP 2 ). Thus, the area can 
always be calculated by summing up over local contribu- 
tions from basic building blocks (pixels, e.g.). Finally, the 
area of a convex body can be continuously approximated 
by the areas of a sequence of convex polygons (conditional 



continuity of Vb). 

There arc other geometric descriptors that share these 
properties with the area. The perimeter is a case in point. 
However, the class of motion-invariant, additive and con- 
ditionally continuous descriptors is not unbounded. Let 
us point this out in full generality for d dimensions. Con- 
sider an arbitrary pattern P that can be decomposed into 
a set union of finitely many convex bodies. According to 
Hadwiger's characterization theorem |27I28| there are only 
(d + 1) linearly independent measures Vq(P), Vd+i(P) 
that obey motion-invariance, additivity and conditional 
continuity. They are called (scalar) Minkowski function- 
als. Thus, in our case of d = 2, the area Vb, the perimeter 
4Vi and the Euler characteristic V2 constitute a complete 
family of scalar morphological measures. Note, that the 
Euler characteristic is a topological invariant and equals 
the number of connected components minus the number of 
holes for patterns in K 2 . The Minkowski functionals were 
applied to neuronal cell classification in [171118) . 
Like the area Vb, the perimeter V% and the Euler char- 
acteristic V2 can be decomposed into local contributions. 
This time they arise from the boundary dP of the body 
P. For smooth boundaries one has 

Vi = \ I d 1 ^, V 2 = ^ / c&S, (2.1) 
4 Jop 27r Jap 

where c denotes the curvature of dP and varies as one 
moves along dP. The factor | is a pure matter of conven- 
tion. For pixel sets, which do not have a smooth boundary, 
V\ and V2 can be calculated by summing up contributions 
from the bonds that confine the pixels, and the corners, 
see P] . 

A natural way of generalizing the concept of the Minkowski 
functionals is to replace the requirement of motion invari- 
ance by motion covariance. Motion covariance means that 
the Minkowski valuations obey simple transformation rules, 
when the body is moved in space: they transform exactly 
as vectors or tensors do under transformations of a coor- 
dinate system. 

The class of motion-covariant, additive and conditionally 
continuous descriptors can be completely characterized by 
a generalization of Hadwiger's theorem |2()U29| . It turns 
out that they can be reconstructed as moments of the 
Minkowski functionals. 

In two dimensions there are three first-order moments of 
the Minkowski functionals, the so-called Minkowski vec- 
tors. For bodies with a smooth boundary, they can be 
represented as follows: 

V = / xd 2 A, Vi = - I xd^ 

JP 4 Jgp 

V 2 = -L / cxd^, (2.2) 
27T Jap 

where x denotes the position vector of the area (perime- 
ter) element d 2 A (d 1 S) to be integrated over. Minkowski 
vectors can also be defined for pixelized images, which lack 
a smooth boundary. 

For the purposes of our analysis, it will be useful to nor- 
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malize the Minkowski vectors and to consider the cen- 
troids: 



Pi = Vi/Vi (* = 0,1,2 if Vt^O). 



(2.3) 



The centroids specify where some aspect of the geometry 
(area, perimeter, curvature) is concentrated. Note, that 
the centroids may, but need not coincide with each 
other. It can be shown that all centroids coincide for spher- 
ically symmetric bodies. 

Moving to second-order moments yields the second-rank 
Minkowski tensors. They are built upon the symmetric 
tensor product denoted by x <S> x =: xx =: x 2 . In two di- 
mensions there are more than three tensors, because, for 
dP-integrals, instead of calculating moments with respect 
to the spatial position x, one may also consider the local 
normal n of the boundary, which points outwards and is 
normalized to one. 1 Thus, for the integrals J^p^S and 
J gp cd*S three types of second-order weights for building 
moments are available, viz. x r n s , where (r,s) = (2,0), 
(1,1) and (0,2) (since we only consider symmetric mo- 
ments, nx and xn are identical). Thus, altogether the fol- 
lowing seven tensors can be formed: 



V<?'°=J xx d 2 A, 
4 JdP 



V r 



1 

27 



op 



(2.4) 
(2.5) 
(2.6) 



In practice, however, we need not consider all of these ten- 
sors, because some of them are linearly dependent |21| . It 
can be shown that only the following tensors carry inde- 
pendent information: 

y 2 ' = f xx d 2 A (2.7) 
Jk 

V?>° = \[ xxd 1 ^ V 1 °' 2 = \[ nnd 1 ^ (2.8) 

4 JdK 4 J OK 

(2.9) 



V4 



I OK 
exxd 1 ^. 



OK 



In the following we will concentrate on these tensors. They 
are listed together with their names in Tabled The nu- 
merics for calculating the Minkowski valuations for pix- 
elized data sets is described in |3*T] . 

Because of motion covariance, the numerical values of the 
second-rank Minkowski tensors depend on the choice of 
the coordinate system. But in many applications, there is 
a natural choice for the origin of the coordinate system. 
For our neuronal cells we will simply take the position of 
the soma as the origin (in other cases it might be useful 
to calculate the second-rank Minkowski tensors V[' s with 
respect to the corresponding centroid Pi for i = 0, 1, 2). 
In order to illustrate very briefly how the Minkowski val- 
uations work for pixelized data sets, let us consider three 

1 First-order moments regarding the normal vectors always 
vanish, as is shown in |3()| 



Symbol 


Formula 


Name 


V 




area 


Po 


j p xd*A/V 


center of mass 


v o 


/ p xxdM 


mass tensor 


Vi 




length of perimeter 


Pi 


LpXd^/Vi 


center of perimeter 


vf' 1 


/ flP xxd l S 


perimeter tensor 




J flp nnd A S 


n-weighted perimeter tensor 


v 2 


f a?S 


Euler characteristic 


P2 


J ffp c^S/V 2 


curvature centroid 




J dp exxd 1 ^ 


curvature tensor 



Table 1. The Minkowski valuations used in this paper. 



simple toy examples (some more examples can be found in 
[2])- They are shown in Figure The red (medium grey) 
filled square, the blue (dark grey) open square and the 
green (light grey) x denote the centroids Po, Pi and p 2 , 
respectively. The tensors are calculated around the center 
of the black square in the middle of the pixel sets as origin. 
The red (medium grey), blue (dark grey) and green (light 
grey) ellipses within the neurons visualize the tensors V^ 2,0 , 
V] 2,0 and V 2 2 '°, respectively. The ellipse for the tensor V®' 2 
is shown at the left-hand side. The equation defining the 

ellipses is always: x = c + a cos(0)e> +sin(0)e < ^, 

where <f> runs from to 27r, e> (e<) is the eigenvector 
corresponding to the larger (smaller) eigenvalue r> (t<) 
of the tensor and c is the center of the soma (except for 
V] 0,2 ; its ellipse is shifted to the edge of the panels). So the 
axis ratios of the ellipses are the ratios of the eigenvalues, 
and the ellipses point into the direction of the eigenvector 
with the larger eigenvalue. The size of the ellipses does not 
carry specific information because of the free scale factor 
a > 0. 

In the top panel of Figure ^ the pixel set displays an ax- 
ial symmetry and is almost point symmetric. Accordingly, 
the centroids are very close to each other; they fan out 
along the symmetry axis. The tensors V 2 ' align perpen- 
dicular to the symmetry axis, because the whole pixel set 
is more elongated along the horizontal axis. The tensor 
ellipses for the mass tensor Vq'° and the perimeter tensor 
V" 2 ' almost coincide, whereas the ellipse corresponding to 
V 2 ' is a bit more elongated. The reason is that the cor- 
ners, which play an important role for the curvature tensor 
V2' are further away from the middle black square, which 



only contributes to Vq'° and 
For the middle panel, the figure has been slightly modified: 
in order to destroy the symmetry, we rearranged one of the 
"arms" . As a consequence, the average pixel is lower down 
than in the first panel, so all centroids move downwards. 
The effect is most prominent for the curvature centroid p 2 , 
because it depends on corners, some of which disappear 
for the rearranged dendrite. Note, furthermore, that the 
centroids span a non-degenerate triangle, a fact that can 
be taken as indicating asymmetry. The lack of symmetry 
is also reflected by the tensor ellipses, which are not par- 
allel any more. Note, furthermore, that the ratios between 
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Fig. 1. Three toy examples to be discussed as an illustra- 
tion. For the centroids the following point styles are used: red 
(medium grey) filled square: po; blue (dark grey) open square: 
pi; green (light grey) x: p2. The ellipses carry information 
about the Minkowski tensors; for more information about the 
construction of the ellipses see the main text. Red (medium 
grey) ellipse: V^'°; blue (dark grey) ellipse: V 2 ' ; green (light 
grey) ellipse: V 2 '° ' ■ Ellipse at the left hand side: Vj ' 2 . 

the bigger and the smaller eigenvalues are larger for the 
second pixel set. The reason is that - due to the "move- 
ment" of the upper right arm - the vertical extension of 
the pixel set shrinks on average, such that the whole body 
is more elongated. 

The bottom panel shows a variation of the body in the 
middle panel, where two holes have been added. This re- 
sults in an Euler characteristic of —1. There is no major 
effect for both po and pi and the related tensors. But for 
P2 a big jump can be observed, and the ellipse for the 
curvature tensor V^ 2 ' is twisted and more elongated. The 



position of P2 can be explained as follows: The hole at 
the right-hand side makes a big negative contribution to 
V2. So, if V2 is calculated around the center of the black 
square, it points to the left hand side. But since the Euler 
characteristic V2 itself is negative, P2 is bounced back to 
the right hand side due to its normalization through V%. 
For the curvature tensor ellipse there is some kind of re- 
pulsion from the right hole, because this hole makes a big 
negative contribution to the tensor; the effect of the other 
hole is much smaller because it is closer to the soma. 
The tensor V^ 2 is shown at the left hand side. It always 
aligns parallel to the grid axis, the reason for this being 
that it crucially depends on normals that can only point 
into four directions for a square lattice. 2 The shape of the 
V®' 2 ellipse can be understood as follows: The eigenval- 
ues of V 1 ' count the number of bonds with horizontal or 
vertical normals, respectively. For all toy examples there 
are more vertical normals, so the tensor is anisotropic. By 
moving from the top to the middle panel, more horizontal 
than vertical normals are destroyed; in this way the tensor 
becomes even more anisotropic. 

Let us conclude this part by adding two comments. First, 
note, that by considering the eigenvalues of a tensor with 
respect to an origin which is given by the body itself, 
motion-invariance is regained. But does this mean that 
we have been returning to the scalar Minkowski function- 
als themselves? The answer is no. Additivity has been lost, 
because forming eigenvalues is not a linear operation, and, 
as a consequence, the eigenvalues of a Minkowski tensor 
cannot be decomposed in the same way as the area is. So 
we have significantly extended the Minkowski framework 
without having given up its conceptual foundations. 
Second, there is a natural extension of our framework to 
three-dimensional neuron data. 



3 The analysis of pixelized neuron data 

Data. We analyze two-dimensional neuron data made a- 
vailable by the courtesy of Prof. Berson [21]. We have 
pixelized maps often neurons. They are assigned different 
types (a, f3, 5, e, 77, t, k, A, 6*, Q. The neuron maps greatly 
differ in terms of scale. Each neuron can be thought of as a 
subset of filled pixels within a square lattice. Not all of the 
neuron pixel sets are connected; some of them consist of 
disconnected parts. This is probably due to an artifact of 
the neuron observations. We will therefore apply a simple 
smoothing. 

Method. For each cell we construct parallel sets with a 
ball of radius r s on a pixel approximation. The parallel 
set P Ts of a body P comprises all points x such that the 
distance between x and P is r s at most. The smoothing is 
illustrated in Figure [3 where the A-neuron is considered. 
In the sequel, the smoothing length will be varied and used 
as a diagnostic parameter. It serves to probe structures at 

2 For an elementary proof, you can start with a single pixel 
and then use additivity. 
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Fig. 2. Two smoothed versions of the A-neuron. Top panel: 
smoothing length: 2 pixels. Bottom panel: 8 pixels. The pic- 
tures are based on data obtained by |25| (their figure no. 5, 
copyright permission by Nature Neuroscience ) . 

different scales. 

For each neuron that has been smoothed with a particu- 
lar smoothing length, we calculate the scalar Minkowski 
functionals, the centroids and the second-rank tensors. For 
the tensors we choose the center of the soma as a natural 
origin. The soma and its center are identified visually, in 
an interactive way. 

Results. We show the neurons with some of the results 
for a smoothing length of one pixel in Figs . 01 — El 3 

3 In the following, one has to be cautious in interpreting the 
green (light grey) ellipses, because for our neuronal cells, the 




Fig. 3. Neurons of type a (top panel), /3 (middle panel) and 
8 (bottom panel). The smoothing length is one pixel. The 
meaning of the points and the ellipses is explained in Fig. 
The small dash in the upper right corner of each panel has 
a length of 20 pixels. The pictures are based on data ob- 
tained by |25| (their figure no. 5, copyright permission by 
Nature Neuroscience). Note, that in all panels of this figure 
as well as of Figs.|3]and|5]the tensor ellipses for V ' and V 1 ' 
almost coincide. 



Let us start with some qualitative observations. First, the 
centroids po through P2 are typically not within the soma. 
Recalling that the centroids are morphological centers, we 
can equivalently say that the soma is quite often eccentric. 
It would be interesting to know whether the eccentricity 
of the soma is characteristic for some types of neurons (for 
this we would have to investigate larger statistical samples 
of neurons). We suspect that the eccentricities depend on 
the function and the local environment of the cells. Fur- 
ther investigations are needed to explore this effect. 
Second, we observe that typically po and pi almost coin- 
cide, whereas P2 may be further away from them. Some- 
thing similar is true about the tensors: The tensor ellipses 
of Vq'° and Vf' often closely resemble each other, whereas 
the ellipse for V^ 2 ' greatly differs. The reason is as follows: 
As our toy examples have shown, p 2 , V^ 2 ' and the corre- 
sponding Minkowski functional (viz. the Euler character- 
tensor V2' sometimes has one or two negative eigenvalues. In 
this case, the ellipse will become smaller and point into the 
direction of e< instead of e>, if |r>| < |r<|. 
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Fig. 4. Neurons of type e (top panel), rj (bottom left panel) 
and l (bottom right panel). The smoothing length is one pixel. 
The pictures are based on data obtained by |25| (their figure 
no. 5, copyright permission by Nature Neuroscience). 



istic) are sensitive to holes. For positive Euler characteris- 
tics, every hole that is off-soma pushes P2 onto the other 
side of the soma. As a consequence, the location of P2 and 
the form of V 2 ' very much depend on the holes, their 
forms and positions. The holes in turn depend on tiny de- 
tails of the branching structure that are not reflected in 
Po and pi and the corresponding tensors V 2,0 and V] 2,0 . 
- Note, that most of the holes are probably due to the 
projection of the neuron into two dimensions. 
We will now turn to a more quantitative analysis. We will 
show several morphological characteristics that are based 
upon the Minkowski valuations as a function of smoothing 
length r s . The point styles designating the different kinds 
of neurons are explained in the top panel of Fig. 
We show the first scalar Minkowski functional Vb for a 
large range of smoothing lengths r s in the bottom panel 
of Figure |H| For very small r s , Vq grows very quickly, as 
r s increases; whereas for larger smoothing lengths, a more 
moderate growth can be seen. For some neurons it ap- 
pears to be linear, for other cell types the function Vo(r s ) 
is clearly convex in this range. Bigger neurons typically 
grow faster than smaller ones. The explanation is as fol- 
lows: Let us consider the (3 cell first. Its overall shape is 
roughly spherical, and its extension 2tq is about 50 pixels. 
If the (3 cell is smoothed with a very large r s > tq = 25, 
all of its substructure is washed out, and we have ap- 
proximately the same result as if a circle of radius tq was 
smoothed by r s . So the volume is about Vq rs Tr(ro+r s ) 2 = 




Fig. 5. Neurons of type k (top panel), £ (middle left panel), 
9 (middle right panel) and A (bottom panel). The smooth- 
ing length is one pixel. The pictures are based on data ob- 
tained by 1251 (their figure no. 5, copyright permission by 
Nature Neuroscience). 



ttTq + 2irror s + irr 2 , which is parabolic in r s . For r s < tq, 
the linear term 2irrQr s is most significant, so the function 
Vo(r s ) appears to be linear in a certain range. 
More generally, let CP denote the convex hull of a pix- 
elized data set P (or, more precisely, the pixel approxi- 
mation of its convex hull). For large smoothing lengths, 
the parallel bodies of P and CP, P rs and CP rs are very 
close to each other; consequently the difference Vb (Pr s ) — 
Vi (CP rs ) is small compared to Vb (Pr B )- The size of the 
parallel body CP ra can be calculated using Steiner's for- 
mula (see p. 367, e.g.): 



V (CP rs ) = V (CP) + r s AVx{CP) + nrj 



(3.10) 



This again defines a parabola, where the Minkowski func- 
tionals Vb and V\ of CP arise as coefficients. As a conse- 
quence, if r s is large enough, the volume Vo(P rs ) is largely 
determined by the Minkowski functionals of the convex 
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Fig. 6. Top panel: the point styles to be used for the differ- 
ent types of neuronal cells. Bottom panel: the volume Vo as a 
function of r s for a large range of revalues. 



hull CP. For small neuronal cells such as the (3 neuron, 
this behavior sets in quite early. Bigger neurons will have 
larger values of Vb(CP) and V\{CP) such that their area 
Vo is larger. 4 

In order to observe the fine-grained structure of the cells 
where the neurons significantly differ from their convex 
hull, we have to concentrate on smaller smoothing lengths 
r s < 20. In Figure [7| the scalar Minkowski functionals are 
plotted vs. the smoothing length r s . For most neurons, 
initially, Vq grows comparatively quickly; around r s = 5, 
however, the growth slows down. As a reason, the arms of 
the neurons that have been blown up, when the parallel 
set was constructed, start to overlap with each other, such 
that increasing r s will not necessarily fill many pixels that 
have not yet been occupied so far. 

For some bigger neurons (a, 8, k, e.g.) a kind of crossover 
can be observed around r s = 5. For the other types of 
neuronal cells, the crossover is less pronounced. 5 
We will now consider V\. For small r s , V\ decreases as a 
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4 Similar considerations apply to V\ and V%. 

5 Note, by the way, that there are plateaus at the zero 
points for the V% vs. r 3 curves. More generally, these curves are 
not continuous, but change stepwise because of our pixelwise 
smoothing. This can be seen, if the r s resolution is enhanced. 
In the following we will neglect discontinuities of this kind; 
they are a pure artifact of our smoothing and do not carry any 
physical meaning. 



Fig. 7. The scalar Minkowski functionals as a function of the 
smoothing length for all cells. (Note, that in the first panel the 
curve for the 6 cell is between the curves for the rj and £ types.) 



function of r s , since V\ is dominated by small scale fea- 
tures that are smoothed away stepwise. V\ reaches a con- 
stant value later on. This is not what one would expect 
for a convex body. The reason, of course, is that the figure 
is far from being convex: As r s increases, V\ will gain at 
the outer parts of the cells, but loose in the inner parts, 
because holes are being filled. Gains and losses roughly 
compensate each other. Note, that the curves for the a, 
S and k cell type show an inflection point, which very 
roughly coincides with the position of their crossover in 
V . 

The curves for the Euler characteristic display a number 
of discontinuities, but there is some more general pattern. 
The negative values indicate that the cells are dominated 
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by holes. For the bigger cell types (a, 5, e and k) there is 
a characteristic dip for small r s . Up to this point, addi- 
tional holes are formed, as branches of the neuron start to 
touch each other. The minimum of the dip roughly seems 
to coincide with the point where Vq shows the crossover 
for the bigger neurons. Similar results for the Minkowski 
functionals have been obtained in |17|. 
A useful way of combining the information present in the 
scalar Minkowski functionals is to construct the following 
dimcnsionlcss quantity Q: 



Q 



ttVo 



(3.11) 



This is a variation of the so-called isoperimetric ratio. For 
a convex body P we have Q{P) > 1, where the equal- 
ity holds for a circle [32,33,34 . Q is considerably larger 
than one, whenever the body under investigation has an 
"excess perimeter" as compared to its area. We show the 
logarithm of Q as a function of r s in Fig. [S] The interpre- 
tation is as follows: For small smoothing lengths r s , most 
of the dendrites are still present; they produce huge excess 
areas for which reason Q starts with very high Q-values. 
As the smoothing length increases, Q goes down. The a, 
e and k cells have the largest Q-values, whereas the /3-cell 
has the lowest Q-values for a large range of smoothing 
lengths because of its smallness and its overall spherical 
shape. For r s < 6 the decrease in log 10 (Q) seems roughly 
to be linear, where the slopes vary with the cell type. 6 
In Figs. R31 and 1101 we consider the centroid distances Pi- 
soma, dis^. For i = they are relatively stable as a function 
of r«, whereas for i = 1 more variation can be observed. 



6 Note, that "linearity" holds only up to discreteness effects 
due to our pixelwise smoothing. 
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Fig. 9. The distance soma - po (top panel) and soma - pi 
(bottom panel) as a function of the smoothing length. 



How is this to be explained? Look at the k neuron as an 
example (Fig. 0. In the lower half of the image the dis- 
tribution of small arms is a bit denser than in the upper 
half. Consequently, for small r s , there is a significant con- 
tribution to the perimeter from this part, and this is also 
reflected in the position of pi, which is the center of the 
perimeter. For larger r s « 10, however, the lower, denser 
part is filled more quickly, whereas in the upper part quite 
big holes are left, which contribute to the perimeter, such 
that the position of pi moves upwards. In this way the 
curve for disi contains very detailed information about 
the morphology of the neuron. 

In terms of diso and disi the soma is most eccentric for 
the e neuron. This is also reflected in our visual impres- 
sions. It might be useful, however, to normalize the dis^ 
parameters by some estimate of the cell size. If we would 
do so, smaller cells would have a reasonable chance of hav- 
ing bigger eccentricities. 

For i = 2 fFig. H0|) we observe even larger variations of the 
centroid distances. Plateaus alternate with jumps that can 
ultimately be traced back to discontinuities of the Euler 
characteristic. For small neurons, such as the (5 type, how- 
ever, there is not much variation, because the cell is very 
small and gets completely filled pretty soon. For the a, 6 
and e-type, there is a common pattern: As the smooth- 
ing length increases, the jumps become larger. The reason 
is probably, that for larger smoothing lengths only a few 
holes will appear far off centered. When one of these outer 
holes disappears, P2 jumps considerably. 
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Fig. 11. The anisotropy parameters derived from the mass 
tensor V 2 ' (top panel) and the perimeter tensor Vj 2,0 (bottom 
panel) as a function of the smoothing length. 



In Figures ITTlthroughll3lwe consider the anisotropy of the 
cells. In order to quantify anisotropy we take the eigenval- 
ues of the tensors Vf' , t> and t< and calculate the quan- 
tity anis := 2(r> — t<)/(|t>| + |r<|) < 2. The anisotropy 
parameters derived from different tensors focus on differ- 
ent kinds of anisotropy (the area elements belonging to 
a body might be distributed differently from those of its 
perimeter elements, for instance). As can be seen from 
Fig- El the anisotropies in V ' and V 1 ' are quite stable; 
most often they decrease slowly, as the smoothing length 
increases. This indicates that the cells display large-scale 
anisotropies that are not destroyed by smoothing the cell. 
For some cells (ry, K, £) the anisotropies are considerable. 
For each cell type the anisotropies of area and perimeter 
do not differ greatly. The V^ 2 ' tensor is a bit more sensi- 
tive to small-scale variations of the morphology, however; 
so the anis(V L 2 '°)-r s curves appear less smooth than the 
anis(V r 2,0 )-r s curves. On the other hand, across the range 
of cell types, the variation is quite high. Thus anisotropies 
seem to have a significant discriminative power. 
It is different with the tensor V^ 2 ' , which is considered in 
Figure IT^l The anisotropy derived from this tensor jumps 
back and forth and sometimes reaches values that exceed 
those derived from the other tensors. This performance 
should not come as a surprise, since we have already seen 
that other characteristics that are related to the Euler 
characteristic such as V^ 2,0 typically show discontinuities. 
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At some point, however, when the smoothing has pro- 
duced one connected pattern without holes (visible for the 
P cell, e.g., where this point is reached quite early), the 
anisotropy stabilizes at a constant value. Apart from this, 
the dependence on r s looks rather chaotic; so far, we are 
not able to extract information that might help to dis- 
criminate between the different cell types. 
As mentioned before, on the square lattice, the last ten- 
sor to be considered, V] ' 2 , has a simple interpretation. It 
checks whether the majority of normals are parallel to the 
horizontal or to the vertical grid axis. If dP is dominated 
by vertical or horizontal normals, V^ ' 2 will display a corre- 
sponding anisotropy; if not, V x ' will roughly be isotropic. 
In Figure lTSl we show some results for single neurons. One 
can learn from this that the anisotropies arising from V^ ' 2 
are quite small. The anisotropy is comparatively large for 
the rj type cell, because this cell is clearly elongated. For 
small values of the smoothing length, anis(y L ' 2 ) is not 
so much influenced by the overall shape of the neuron, 
but rather by the directions of the single arms. Interest- 
ingly, the graphs shown are qualitatively different for the 
different types of cells: One cell (viz. the a cell) starts 
with zero anisotropy, whereas others begin with a non- 
zero anisotropy. Moreover, there are significant peak struc- 
tures. But because of its relation to normals, n, the value 
of V] ' 2 depends to a large extent on the orientation of the 
cell with respect to the grid. For this reason V] ' 2 is only 
of limited use. 



In Fig. El the traces of the tensors V 2,a are considered 
(the trace of the fourth tensor, V^ 0,2 need not to be taken 
into account at this point, because it equals Vi). Quali- 
tatively, the viewgraphs for V 2,0 resemble the curves of 
their scalar counterparts, Vi for i = 0, .., 2. In order to ex- 
tract more specific information, it is thus useful to divide 

Tr (y 2 '°^ by Vi, respectively, for i = 0, .., 2. The result is 
a measure of how concentrated a cell is in terms of area, 
perimeter or curvature: Tt(Vq 2 ' ) /Vo, for instance will be 
the bigger, the further the soma and those parts of the cell 
that bear most of its volume lie apart. Results can be seen 
from FigureEH ^((Vq' )/^ increases continuously, as r s 
is enhanced. The reason is that more and more pixels are 
added at the outer parts of the neuron, so the neuron be- 
comes less and less concentrated. In T!r((Vi ' )/Vi there is 
a kink at least for some neurons (a, 5, e, n). It indicates an 
additional growth effect. Very probably the explanation is 
that for small r s , the small branches within the cell signifi- 
cantly contribute to V^ ' , so the neuron appears to be very 
concentrated; for larger values of r s the arms merge and 
do not contribute to the perimeter any more, so most of 
the neuron's perimeter is found at its outer parts. Note, 
that the kinks roughly set in at the relocations of the 
crossover point in Vo and the inflection point in V\ for the 
a, 5 and n cells. 



Constructing global measures A multiscale analysis like 
that presented in this paper, leads to rich and detailed 
information on the geometrical aspects of an object. Nev- 
ertheless, once such a description of the data has been ob- 
tained, it is often useful to derive a compact set of global 
measures that summarize the most important morpholog- 
ical aspects. In this paper, we consider several ways of 
condensing multiscale information, i.e. a function of some 
scale, into simple parameters: The monotonicity index [171 
IT%] is defined as 

(3.12) 



s + d + p ' 

where s, d and p count each time the function increases, 
decreases and remains unchanged, respectively. Thus i s 
quantifies the fraction of the interval where the function 
is monotonically increasing. The mean value is the aver- 
age value of the function over the interval. The half scale 
is the scale at which the area below a curve reaches half 
of its total value. A different way of constructing global 
parameters is to consider the slope of some characteristic 
in some particular range of revalues. 
In Figure ITrll we visualize the average slopes of Vo in the 
range r s £ [10,60] and of log 10 (Q) in the range r s £ 
[1.5,4.5] for the different cells. In both cases we choose 
a range of r s -values for which the functionals under inves- 
tigation look roughly linear for most cell types. Results 
are shown in Figure [TBI One can immediately see that the 
slopes discriminate amongst the different cell types. 
In order to further illustrate our approach, we selected 
two two-dimensional feature spaces, which are spanned by 
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Fig. 14. The traces of the tensors ' for i 
function of the smoothing length. 



size- independent morphological characteristics. In order to 
calculate them, we considered the interval r s <E [0, 20] and 
a spacing of 0.2. 

Our first feature space is spanned by the mean of the 
anisotropy parameter derived from Vq'° , anis ^V^ 2 ' ^, and 
the mean of the anisotropy parameter corresponding to 
Vq'°, anis ^Vf' ^ . It is shown in Figure nTf a). There ap- 
pears to be some systematic correlation between both 



Fig. 

cells 



13. The anisotropy parameter derived from Vf' 2 for four characteristics: cells with higher anis (V^ 2,0 ) tend to have 

~* ri *-i Inn/tf inn + h n l~i v-v-\ nnf Ininrr I i-r-f- V ' 



as a function of the smoothing length 



higher anis ^V^ 2 ' ^ as well. Given the meaning of these 

characteristics, this should not come as a surprise, al- 
though it is in principle possible to have high anisotropy 
in Vq'° and low anisotropy in V^ 2 ' . Thus, for discriminat- 
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po relatively to their soma. Figure Hlsl presents a dendro- 
gram obtained by a simple hierarchical agglomerative clus- 
tering j35j of the scatter plot distribution shown in Fig- 
ure Etb)- Such a structure suggests a possible taxonomy 
for the ten types of cells. As expected, the cells (3, 77, 9 and 
C are similar, inhabiting the same branch at the lower part 
of the dendrogram. For the remaining subset, the cells a 
and k end up markedly distinct from the group of cells 
formed by 5, e, 1 and A. 

Although the proposed methodology may have a bearing 
on the classification of cat ganglion cells, it is difficult to 
make more definitive conclusions at this point, because 
the original classification 125] takes into account not only 
the neuronal morphology, but also the cell stratification 
and the size of the soma. Moreover, except for the more 
common a and j3 types, only a small number of examples 
of the cell types have been analyzed in the related liter- 
ature |26j. A more detailed examination of which feature 
spaces are most useful has to wait for further data. 



4 Conclusions 



ing between different cells, one dimension of this feature 
space is essentially redundant. But the presence of some 
correlation might be used to describe some common trait 
shared by all cells. 

A different situation can be observed for our second fea- 
ture space. It is spanned by the monotonicity index i s (diso) 

and by the half scale h (Tr (^i 2 '°)) ■ As can be seen f rom 
Figure ITTT b). the scatter is larger, and the cells form kind 
of groups. Note, in particular, that neuronal cells that look 
similar at least in some respect tend to appear close to 
each other in this scatter plot. For instance, cells j3, i], 8 
and £ are close to each other in the bottom panel of Fig- 
ure El especially regarding the position of their center of 
mass 



We have analyzed two-dimensional projections of neuronal 
cells using higher-order Minkowski valuations. Our mea- 
sures detect different kinds of substructures, providing a 
natural extension of previous works that deal with the 
more traditional shape functionals |17II18| . An extensive 
discussion of the results obtained for a set of ten neu- 
ronal cells was included that illustrates the interpretation 
of the suggested measures and implications for neuromor- 
phometric studies. As far as our limited set of samples is 
concerned, significant similarities and differences between 
the cell types have been found, leading to a putative tax- 
onomy. It is a pending question whether the differences 
found will still be characteristic of the types in a statisti- 
cal sense. 
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Fig. 17. Scatter plots from selected features of an extended 
Minkowski analysis showing the population of the feature space 
with the neuronal cells. 
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